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ABSTRACT 

The ellipticities and alignments of clusters of galaxies, and their evolution 
with redshift, are examined in the context of a A-dominated cold dark matter 
cosmology. We use a large-scale, high-resolution N-body simulation to model 
the matter distribution in a light cone containing ~ 10^ clusters of mass M > 
2 X 10^^ Mq out to redshifts oi z = 3. The best-fit three-dimensional ellipsoid 
of the mass distribution is determined for each cluster, and the results are used 
to analyze cluster ellipticities as a function of mass, radius, and redshift. A 
similar analysis is done in two dimensions in order to allow direct comparisons 
with future observations. Cluster ellipticities are determined within different 
radii, including 0.5, 1.0, and 1.5 comoving Mpc. We find strong cluster 
ellipticities: (e) = (1 — a2/ai) ~ 0.3 — 0.5. The mean ellipticity increases with 
redshift from (e) ~ 0.3 at 2; = to (e) ~ 0.5 at z = 3, for both 3D and 
2D ellipticities; the evolution is well-fit by (e) = 0.33 + 0.05z. The ellipticities 
increase with cluster mass and with cluster radius; the main cluster body is more 
elliptical than the cluster cores, but the increase of ellipticities with redshift 
is preserved. Using the fitted cluster ellipsoids, we determine the alignment of 
clusters as a function of their separation. We find strong alignment of clusters for 
separations < 100 Mpc; the alignment increases with decreasing separation 
and with increasing redshift. The evolution of clusters from highly aligned and 
elongated systems at early times to lower alignment and elongation at present 
refiects the hierarchical and filamentary nature of structure formation. These 
measures of cluster ellipticity and alignment will provide a new test of the current 
cosmological model when compared with upcoming cluster surveys. 

Subject headings: cosmology: theory — galaxies: clusters: general — large-scale 
structure of universe 
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1. Introduction 

In this paper we examine two aspects of cluster structure that have received some, 
but not extensive, attention. The first is the distribution of cluster ellipticities and the 
evolution of this distribution with mass and redshift. The second is cluster alignment: how 
aligned are clusters with their neighbors, and how does the alignment change with cluster 
mass, separation, and redshift? Ongoing optical and X-ray surveys, and surveys based 
on weak gravitational lensing and the S-Z effect, will provide larger samples of clusters 
than currently available; these surveys will enable the measurement of cluster ellipticities 
and alignments. When compared with theoretical expectations as provided below, these 
measures, and their evolution with redshift, will provide a test of the cosmological model 
and help our understanding of cluster formation and evolution. 

Observations of galaxy clusters suggest that their shapes, as traced by the distribution 
of galaxies or X-ray emission within the cluster, are not usually spherical but instead tend to 
be elongated (Carter & Metcalfe 1980; West 1989; Plionis, Barrow, & Frenk 1991; de Theije, 
Katgert, & van Kampen 1995; Basilakos, Plionis, & Maddox 2000). Melott, Chambers, 
& Miller (2001) and Plionis (2002) studied clusters to z < 0.13 and suggest evolution in 
the gross morphology of clusters in this range, but Rahman et al. (2004) have pointed out 
that different cluster samples do not give consistent rates of evolution. Floor et al. (2003) 
investigated simulated clusters at 2; < 0.2 and found little evolution in this narrow range; 
Jing & Suto (2002), Suwa et al. (2003), and Rahman et al. (2004) found evidence of ellipticity 
evolution in ACDM simulations out to z = 0.5 — 1. The evolution of cluster ellipticities to 
higher redshifts has not been investigated. 

Conventional models suggest that galaxy clusters form through a process of hierarchical 
clustering, where gas and galaxies flow into denser regions along interconnecting large-scale 
filamentary structures (e.g., Shandarin & Klypin 1984; Gott, Cen, & Ostriker 1996; Ostriker 
& Cen 1996, and references therein). As a result of these inflows, we expect clusters to be 
aligned with their neighbors, especially when both are members of the same filament. Pre- 
vious observations of cluster alignments have yielded controversial results. Binggeli (1982) 
found that Abell clusters tend to ahgn with their nearest neighbor out to separations of 
~ 15 Mpc, and that cluster orientation is related to the distribution of all surround- 
ing clusters out to separations of ~ 50 h"^ Mpc. However, contradictory results have been 
claimed by others; for example, Struble & Peebles (1985) extended Binggeh's method and 
failed to find significant alignment. Fhn (1987), Rhee & Katgert (1987), West (1989), and 
Plionis (1994) found optical alignment for scales below 30/i~^Mpc. X-ray emission from 
clusters has also been used to determine cluster shapes and alignments. Ulmer, McMillan, 
& Kowalski (1989) studied 46 Einstein clusters and did not find alignment; but Chambers, 
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Melott, & Miller (2000) re-examined their results using updated cluster positions and found 
alignment among near neighbors. Rhee, van Haarlem, & Katgert (1992) combined X-ray 
and optical data, and found that members of the same supercluster tend to be aligned. 
West, Jones, & Forman (1995) found correlation among cluster X-ray orientations out to 
~ 10 h-^ Mpc, and Chambers, Melott, & Miller (2002) used Emstem and ROSAT data to 
find alignment for separations below 30 Mpc. 

Advances in computational power enable large-volume, high- resolution cosmological 
simulations that provide a large sample of clusters for studying ellipticities and alignments 
and their evolution with rcdshift. Previous simulations by Splinter et al. (1997), Faltcnbacher 
et al. (2002), and Onuora & Thomas (2000) used 128^256^ and 512^ dark-matter parti- 
cles, respectively, to study elhpticities and orientations of nearby {z ~ 0) clusters; all found 
significant alignments for nearby clusters. In this paper we use a large-scale high-resolution 
simulation of the current best-fit ACDM cosmological model, with 2 x 10^ particles, to 
study cluster ellipticities, orientations, and alignments in much greater detail than previ- 
ously possible. Wc examine cluster alignments as a function of cluster pair separation and 
find significant structure in this relation for separations < 100 Mpc, and also determine 
the evolution of cluster ellipticity and alignment with redshift, from 2; = to 2; = 3. This 
work seeks to provide a foundation for comparison with future observations of large cluster 
samples to high redshifts. Such comparison can be used to test the predictions of the current 
cosmological model as well as enhance our understanding of cluster formation and evolution. 
In order to provide easy comparison with observations, all important quantities, such as 
ellipticity and alignment measures, are obtained using both the three-dimensional halo dark 
matter particle distribution and the projected two-dimensional particle distribution; the lat- 
ter provides a direct comparison with observations. To further facilitate comparison with 
observations, our analyses of ellipticity and cluster alignment are repeated using fixed cluster 
radii typically used in observations, rather than the friends-of-friends (or virial) radius. 

We describe the simulation and cluster selection in §2. The derived cluster ellipticity 
distribution and its dependence on cluster mass and redshift are presented in §3. The 
alignment of clusters as a function of cluster separation and mass, and the evolution of 
alignment with redshift, are discussed in §4. Our conclusions are summarized in §5. 

2. The Simulation 

Wc use a high-resolution N-body simulation of a fiat universe with matter density flm = 
0.27 (baryon density Q\^h^ = 0.023 plus cold dark matter), cosmological constant = 0.73, 
Hubble constant Hq = /i lOOkms"^ Mpc~^ = 70kms~^ Mpc~^, primordial scalar spectral 



-4- 



index Ug — 0.96, and linear matter power spectrum amplitude — 0.84. These values are 
consistent with both the cosmic microwave background fluctuation spectrum of WMAP as 
well as with large-scale galaxy and cluster results (e.g., Bahcall et al. 1999; Spergel et al. 
2003; Tegmark et al. 2004). 

The simulation contains 1260^ ^ 2 biUion particles in a periodic cube. The box length 
is 1500/i~^Mpc and the individual particle mass is 1.264 x lO^^/i~^M0. The spline kernel 
softening length is e = 17/i^^kpc; thus, the spatial resolution is comfortably smaller than 
the ~ lOO/i^^kpc characteristic core size of clusters, and the high mass resolution ensures 
that two-body relaxation is unimportant in the cluster cores. The initial conditions were 
generated using the publicly available codes GRAFIC2 and LINGERS^ (Bertschinger 2001; 
Ma & Bertschinger 1995). The simulation was carried out using the TPM (Tree-Particle- 
Mesh) code^ (Bode & Ostriker 2003), using a PM mesh of size 1260^. The opening angle 
in the tree portion of the code was 9 — 0.577 ^ 3~^/^. The time step parameter was set 
to ?7 = 0.3, and the initial domain decomposition parameters were A — 2.0 and B = 12.0 
(for details on these see Bode & Ostriker 2003). The computation was carried out using 420 
processors on the Terascale Computing System at the Pittsburgh Supercomputing Center, 
and the running time was approximately five days. At the end of the simulation ~ 10^ 
trees — occupying 2% of the volume (including all PM cells with an overdensity above 39) 
and containing 40% of the mass — were being followed at full resolution. In the following, 
only objects containing 160 or more particles, i.e. with masses above 2 x 10^^ h^^ Mq, will 
be considered; this is four times above the hmit imposed by the densest PM-only cell, thus 
we can be confident that resolution will not cause a problem in our analysis. This mass 
threshold corresponds to typical groups of galaxies, and thus includes all clusters, from poor 
to the richest systems. 

Particles were saved at each PM time step in order to later reconstruct the mass distri- 
bution along the light cone. For distances less than 1500 Mpc from the simulation origin 
at ^ = 0, all particles in spherical shells defined by each PM time step were saved (although 
only for distances less than half the box size, or 2; < 0.265, is the shell contained within the 
simulation volume). For distances above 1500 hr^ Mpc, a single octant with all three spatial 
coordinates greater than zero is saved. Data were saved for a total of 670 redshift shells in 
the light cone between z — 3 and z — 0. Since the fight cone is larger than the simulation 
box size, a given cluster can appear more than once, but at different redshifts. Cluster halos 
were identified using a standard FOF (friends-of-friends) halo finding routine, run on the 
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light cone output, using a linking length parameter b — 0.2 times the mean interparticle 

separation, corresponding to an overdensity of 180 times the mean at all redshifts (Lacey & 
Cole 1994; Jenkins et al. 2001). Different linking lengths (6=0.16 and 0.25, corresponding to 
fixed overdensities of 360 and 92) were also tried for comparison, yielding consistent results. 
For comparison, the virial overdensity for spherical collapse ranges from 97 times the mean 
at z=0 to 174 at z=3 for the model used in this paper (Bryan & Norman 1998). 



3. Cluster EUipticities and their Evolution 

We determine the structure of halos in the simulation by fitting each identified halo to 
a three-dimensional ellipsoid and recording the lengths and orientations of the three axes of 
the fitted ellipsoid. First, the center of each halo was taken to be the halo center of mass. 
Once the center of mass of a given halo was identified, the best-fit ellipse was found by taking 
the matrix of second moments of particle positions about the center of mass, defined as 

where the sum is over all particles in the halo. The dark matter particles are assumed to trace 
the distribution of galaxies within the cluster, as suggested by observations. All particles 
in our simulation have the same mass, so m was set to 1. The normalized eigenvectors of 
/ correspond to the unit vectors of the axes of the best-fit ellipse, and the corresponding 
eigenvalues (Aj) give the lengths of these axes (oj) by = v^SA^; the factor of v^S is merely a 
geometrical normalization factor, which will cancel out in our analysis. The eigenvalues and 
corresponding eigenvectors were sorted such that Ai > A2 > A3. We use the most common 
measure of cluster ellipticity, defined as 

6 = 1- 02/01 = 1 - VA2/A1, (2) 

where ai and 02 are the primary and secondary ellipsoid axes lengths, respectively. We note 
that other definitions of eUipticity, e.g. e = 1 — (02/01)^ (de Theije, Katgert, & van Kampen 
1995; Splinter et al. 1997) and eccentricity, e = ^1 — (02/01)^, are used occasionally in some 
studies. We highlight these definitions in order to clarify comparisons but, unless otherwise 
specified, we use relation (2) for the definition of e in our analysis, as it has the most obvious 
relationship between ellipticity and axis ratio. Clearly, e and e increase monotonically with 
one another inside the physical range < e < 1, and our results are not qualitatively changed 
by using either definition. We used Monte Carlo simulations to test the algorithm's ability 
to recover halo eUipticities, and find that the algorithm recovers ellipticities to within ±0.01, 
with accuracy increasing with higher ellipticity. The results are insensitive to the number of 
halo particles (with a slight increase in dispersion for a lower number of particles) . 
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In order to facilitate comparison with future observations, the process was repeated by 
viewing halos in two-dimensional projection, as in observations. Because halo orientations 
are, on average, random, we use the projection in the x-y axis. The eigenvalues and corre- 
sponding eigenvectors were again sorted such that X2d,i > ^2(1,2, and the projected cllipticity 
of each halo was taken to be = 1 — 02^,2/02^,1- Figure 1 compares the two and three- 
dimensional ellipticities found for the simulated halos. The two ellipticities are correlated as 
e ~ ^2d] the correlation is stronger with increasing ellipticity. 

We find that clusters have moderate to large ellipticities: (e) > 0.3. The distribution of 
ellipticities peaks in the range 0.3 ^ e < 0.5, corresponding to an axis ratio of | < 0,2/0-1 ^ |. 
Thus clusters are not spherical, but elongated ellipsoids. This holds for clusters of all masses. 
Figure 2 shows the ellipticity as a function of cluster mass in six redshift bins of width 
Az — 0.5. A weak correlation of ellipticity with cluster mass is observed; high-mass clusters 
tend to have high eUipticities, but the scatter is large. This is investigated in more detail 
below. 

The evolution of cluster ellipticities with redshift is illustrated in Figure 3. The results 
for the 3-D and 2-D eUipticities are presented by the solid and dashed histograms, respec- 
tively. The dotted vertical line shows the mean 3-D ellipticity at the given redshift. The 

results reveal that the mean ellipticity increases with redshift from (e) ~ 0.3 at low redshift 
to (e) ~ 0.5 at 2; ~ 2.5. It is also clear that despite the dispersion in the €24 — e relation, the 
distributions of 3-D and 2-D ellipticities arc almost identical. This is a consequence of sec- 
ondary ellipticities Cgcc = 1 — 03/02 having a similar distribution to primary ellipticities in this 
range (well approximated by a Gaussian with (cscc) = 0.2 and a^^^^ = 0.12); thus the primary 
ellipticity e is reproduced on average in projection by triaxial clusters with the above mean 
axis ratios. Much of the cluster population has moderate to large ellipticity (0.2 < e < 0.6), 
with larger ellipticities at higher redshift. Although the dispersion of individual ellipticities 
about the mean value is large. Splinter et al. (1997) and Floor et al. (2003) found that, when 
multiple simulations of the same cosmology are run, the mean ellipticity at a given redshift 
has a very low dispersion, (Tmcan ^ 0.03. The best-fit Gaussian to the observed distribution 
of z ^ cluster ellipticities ((e) ~ 0.34, 0"^ ~ 0.16) from Plionis, Barrow, & Frenk (1991) is 
shown for comparison in Figure 3; it is consistent with the 2; ~ ACDM cluster ellipticity 
distribution (with a slight offset due to the larger mass threshold in the observations) . 

The evolution of ellipticity with redshift is illustrated more clearly in Figure 4. The mean 
ellipticity of all clusters as a function of redshift is plotted for the 3-D ellipticities (bold line). 
The mean 2-D ellipticity is not shown but is identical to the 3-D at all redshifts to within 
±0.05. The mean ellipticity rises consistently from (e) ~ 0.33 at 2; = to (e) ~ 0.5 at 2; = 3. 
A least-squares linear fit yields (e) = 0.3276(±0.0003) + 0.0513(±0.0003)^ (rms = 0.0048). 
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For the projected 2-D ellipticities we find (ea^) = 0.3317(±0.0003) + 0.0490(±0.0003)^ (rms 
= 0.0048). Thus we find an approximately linear increase in mean ellipticity from ^ = to 
z = 3 with (e) = 0.33 + 0.05z, for both the 3-D and 2-D ellipticities. This slope is consistent 
with those obtained in previous ACDM simulations at z < 0.5 (Rahman et al. 2004) as 
well as with some observations aX z < 0.1 (Melott, Chambers, & Miller 2001), although, as 
pointed out by Rahman et al. (2004), different observational catalogs disagree on the slope. 

The evolution of ellipticity as a function of cluster mass is also presented in Figure 4. 
Clusters are divided into three mass bins: low-mass clusters with 2 x 10"*^^ < M < 5 x 
10"^^ Mq, medium-mass chistcrs with 5 x 10^'^ < M < 10"*^^ M©, and high-mass clusters 
with M > 10^'^h~^MQ (where cluster mass is the total FOF mass). We find, as expected 
from Figure 2, that higher-mass clusters are more elliptical at all redshifts. The ellipticity- 
redshift relation shows the same trend for clusters of all masses, with higher mass clusters 
shifted to larger mean eUipticity at any given redshift. We also find that the projected mean 
ellipticity traces the mean three-dimensional ellipticity in all three mass bins. The mean 
ellipticity-redshift relation found for all masses is dominated by the low-mass clusters, as 
expected, especially at high redshift where higher-mass clusters become rare. 

The rate of change in eUipticity could be affected by numerical effects in the cores of 
the simulated clusters, particularly when the number of particles in the cluster is small. To 
test if numerical relaxation is occurring, we repeated the analysis with a higher mass and 
force resolution simulation, evolved with the same code. The cosmological parameters for 
this simulation arc fairly similar to the light cone run (Q^n = 0.3, fl^ = 0.7, h=0.7, ns=l, 
and (78=0.95), but the box size is smaller (L = 320/i~^Mpc) and hence the resolution is 
higher: N — 1024^, making the particle mass 2.54 x 10^ Mq; the spline softening length 
is 3.2 kpc. The FOF (6=0.2) finder was run on the entire simulation cube at different 
redshifts to find clusters with masses greater than 2 x 10^^ Mq, (i.e., with at least 7960 
particles per cluster, as compared to 160 in the primary simulation). We find no change 
in our results for high or medium-mass clusters, but do find that the slope for low-mass 
clusters steepens slightly, to ^ ~ 0.06 — 0.07. This leaves our ^ ^ 1 results unchanged but 
implies that ellipticities of low-mass clusters at high redshifts may be larger than those in 
our primary simulation by ~ 5 — 8%. 

The clusters studied so far include all dark matter particles within the friends-of-friends 
(FOF) radius. However, this does not correspond to an easily observed quantity, so we repeat 
the analysis using typical observed radii of 0.5, 1.0, and 1.5 Z;,^^ Mpc (comoving), denoted 
by -Ro.5, Ri.o, and Ri,5 respectively. For these analyses, we use the previously determined 
center of mass (determined from all halo particles) as the center of each halo, and then 
drop all particles further than the given radius of interest. The analysis then follows the 
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same methods as above, but using only the particles within the revised radius. In order to 
ensure that the friends-of-friends halo radius -Rfof extends to the cut-off radius of interest, 
the proper mass threshold must be used for each radius. The initial mass cut, M > 2 x 
10^^ h^^ Mq, is sufficient to ensure that all halos have i?FOF > -Ro.5- Similarly, a mass cut of 
M > 4 X 10^3 h-^ Mq ensures i?FOF > Ri.o, and M > 10^^ h'^ M© ensures i^poF > ^i.5- 

The results are presented in Figure 5. The figure shows the evolution of the mean ellip- 
ticity with rcdshift for both the 3-D and 2-D cases, for cluster radii of -Ro.5! -Ri.Oi and -R1.5. 
We find that the measured ellipticities (both 3-D and 2-D) decrease strongly for smaller 
radii (for the same Mfof clusters). As the halo centers arc the most dense, one expects an 
increasingly spherically symmetric distribution closer to the center, due to the shorter relax- 
ation time. Additionally, extended structures and filaments which would strongly contribute 
to measured ellipticity are increasingly reduced with decreasing radius. At i?o.5, most ellip- 
ticities have decreased by ~ 0.1 — 0.2 from their it!i.5 or -Rfof values. A similar correlation 
between two and three dimensional ellipticities is observed for all radii. The best linear fits 
to the evolution of e with redshift for each cut-off radius yield a slope of ^ ~ 0.05, with a 
somewhat shallower slope of ^ ~ 0.025 for i?o.5- 

To test the dependence of these results on the FOF hnking length b used in the clus- 
ter identification, we identified clusters using two additional different linking lengths and 
repeated the analysis. We find a dependence on linking length similar to the dependence on 
cluster radius, as expected. Using a linking length h = 0.25 (enclosing an overdensity of 92) 
instead of 6 = 0.20 (enclosing an overdensity of 180) yields clusters with larger FOF radii, 
well beyond the virial radius, and thus systematically increases the mean cluster ellipticities 
by ~ 0.1. A linking length of 6 = 0.16 has, as expected, an opposite effect; it systemat- 
ically decreases mean ellipticities by ~ 0.025, similar to reducing the radius to Ri.q. In 
both cases the dependence on redshift is preserved. We note that the fixed radii used above 
do not correspond to the virial scales of clusters, but would be less difficult to implement 
observationally. 



4. Cluster Alignment and its Evolution 

Using the parameters measured in § 3 for each halo, we identify pairs of halos and record 
the separation and alignment measure for each pair. We examine two common measures of 
alignment: first, the angle between the two major axes of a pair of clusters (the "correlation" 
angle), and second, the angle between a given cluster major axis and the line connecting the 
cluster to its paired cluster (the "pointing" angle). We use as our "correlation" alignment 
measure the squared dot product of the two halo major axes (given as unit vectors by the 
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eigenvectors found in § 3), yielding cos^(^c) for each pair, where 9c is the angle between 
the two major axes. We consider cos^(^c) instead of cos(^c) because ±cos(^c) are physically 
identical, since the unit vectors defining the halo axes do not have a preferred "positive" 
direction. We similarly use cos^(^^p) as our "pointing" angle measure, but in this case 6p 
is the angle between a cluster major axis and the line connecting its center of mass to the 
center of mass of its paired cluster. We find, as Splinter et al. (1997), that both measures 
give similar results, with a stronger signal for the pointing angle measure. We also find that 
the pointing angle signal extends to larger scales. 

Using the correlation angle measure, we find that pairs of clusters are strongly aligned 
for separations < 30 /i^^ Mpc. On larger scales, alignment between pairs of clusters is ran- 
dom and uniformly distributed. The alignment function - the mean cluster alignment as a 
function of cluster pair separation r - is presented for different redshift bins in Figure 6 (for 
3-D) and Figure 7 (for 2-D). Assuming a random uniform distribution in 9 {P{9) — |^), one 
expects (cos^(^)) = | for the 3-D case and (cos^(^)) = ^ for the 2-D case (for both correla- 
tion and pointing angle measures); obtaining a (cos^(^)) greater than these values indicates 
alignment. The data yield (cos^(6'c)) = (cos^(^^))ran (where (cos^(6'))ran = \i \ foi^ the 3-D 
and 2-D cases, respectively) for large separations (r > 30/i~^Mpc). At smaller separations 
the data show clear alignment of clusters, with a pronounced rise in (cos^(6'c)) with decreas- 
ing r, to a maximum at r 3 /i"-*- Mpc. The fall in (cos^(^^c)) below ?>h~^Mpc is due to 
the close proximity of the cluster pairs; clusters at such small separations begin to overlap, 
and measures of ellipticity and alignment become less meaningful. The alignment increases 
significantly with increasing redshift. For example, (cos^(^c)) at r 3/i~^Mpc rises from 
~ 0.40 at < ^ < 0.5 to ~ 0.46 at 2.0 < z < 2.5. The signal-to-noise ratio decreases at high 
redshift, as the number of massive clusters decreases. At the highest redshifts, 2.5 < z < 3.0, 
the rise in (cos^(^^c)) with decreasing r (to (cos^(6'c)) ~ 0.48 at r ~ 3/i~^ Mpc) is considerably 
steeper than for z < 2.5. 

Similar trends are observed for the 2-D case in Figure 7. We see (cos^(^c)) — (cos^(^))i.an — 
I for r > 20 h^^ Mpc, with a similar rise in (cos^(^^c)) with decreasing r. Alignment increases 
with increasing redshift, from (cos^(^^c))max ~ 0.54 for < 2; < 0.5 to (cos^(^^))max ~ 0.57 for 
2.0 < z < 2.5. At high redshift, uncertainties are large, and for 2.5 < z < 3.0 the alignment 
function becomes more erratic. In this redshift range, (cos^(^c))max ~ 0.65 at r 2 Mpc. 

Figures 8 and 9 are similar to Figures 6 and 7, but use the pointing angle alignment 
measure described previously. This measure is the angle between a given cluster major axis 
and the line connecting the cluster center of mass to its paired cluster center of mass at 
the given separation. We again use (cos^(^p)) as our alignment measure, employing this 
alternative definition of 9p. We find that the pointing angle alignment function is similar 
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to the correlation angle alignment function, with larger amplitude and a positive alignment 
signal extending to larger scales of ~ 100 Mpc. Previous observations by West (1989) 
and Plionis (1994), which used the pointing angle measure, are shown for comparison in 
Figure 9. The results are consistent within la with the ACDM predicted 2-D alignment at 
z = Q. 

As in § 3, we compare our results with those obtained using the more readily observed 
fixed co-moving radii of 0.5, 1.0, and 1.5 h^^ Mpc. We find that the alignment measures are 
nearly independent (within the la error bars) of the cluster radii used for determining cluster 
orientations. A similar near independence holds when varying the halo linking lengths; 
clusters identified using FOF linking lengths of 6 = 0.25 and 6 = 0.16 yield essentially the 
same results as those using b — 0.20 (with a slight ~ 0.01 increase in (cos^(6')) at small r for 
6 = 0.25). 

We investigate the sensitivity of the alignment function to the masses and elliptici- 
ties of the halos by examining a mass and ellipticity- weighted average of cos^(^), taken as 
{MiMj 008^(6)) /{MiMj) and {e^ej cos^ [9)) / {tit j) , where Mj and Mj represent the masses of 
the two halos in a given pair (likewise for and tj). We find that the ellipticity- weighted 
average closely traces the standard expectation value of cos^(^) at all redshifts for both 
the 3-D and 2-D cases, and thus do not further explore a possible ellipticity dependence of 
this relationship. We note that Onuora & Thomas (2000) also found that highly elongated 
clusters show no greater tendency to alignment than more spherical clusters. 

We do find some deviations between {MiMj cos"^ (9)) / (MiMj) and (cos^(^)), and thus 
explore the possible mass dependence of this relation by dividing clusters into three mass 
bins using their FOF masses: low-mass clusters with 2 x 10^^ < M < 5 x 10^^ H^^Mq, 
medium- mass clusters with 5 x 10^^ < M < W^H-^Mq, and high -mass clusters with M > 
10^^ h^^ Mq. Because alignment is determined using pairs of clusters, each cluster is labeled 
with a letter corresponding to its mass: H (high-mass cluster, as defined above), M (medium- 
mass), or L (low-mass). Pairs are then identified as either H-H, H-M, H-L, M-M, M-L, or 
L-L, with the letters corresponding to the two elements of the cluster pair. We repeat our 
alignment analysis as a function of redshift for each of these six cluster subsamples, using 
both the correlation angle and pointing angle measures of alignment. 

Figure 10 presents the results of this alignment analysis as a function of pair separation 
and redshift for the three equal-mass pair subsamples, using the pointing angle measure. The 
results show that the cluster ahgnment function increases and shifts to larger scales as the 
mass of the clusters increases from L-L to M-M to H-H. The scale at peak alignment of the 
L-L, M-M, and H-H pairs increases from r ~ 2 Mpc to ~ 3 Mpc to ~ 6 h^^ Mpc, re- 
spectively. This is unsurprising, as it suggests that more massive cluster pairs show stronger 
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alignments to larger scales. Mixed mass pairs, H-M, M-L, and H-L, follow intermediate rela- 
tionships, with smaller peak alignments. It appears that clusters of comparable masses yield 
the strongest alignment; the alignment strength and scale increases with increasing cluster 
mass. Results for the 2-D cluster alignments arc presented in Figure 11. The correlation 
angle results are qualitatively similar, but with weaker amplitude and no significant signal 
for separations above ~ 50 Mpc. 

The cluster alignment increases for all masses with increasing redshift; high redshift 
clusters are more strongly aligned than low redshift clusters. The alignment increases more 
rapidly at high redshift among high-mass clusters than among low-mass clusters. For 1.0 < 
z < 1.5, we find a maximum value of (cos^(6'p)) ~ 0.56 for L-L pairs, and a maximum 
value of (cos^(6'p)) ~ 0.65 for H-H pairs. At 1.5 < 2; < 2.0, the maximum alignment of L-L 
pairs rises to (cos^(^p)) ~ 0.58, and the maximum ahgnment of M-M and H-M pairs rises 
to (cos^(^p)) ~ 0.62 (there are too few H-H pairs for statistical analysis). At 2.0 < 2; < 2.5, 
the L-L pairs have a maximum of (cos^(^p)) ~ 0.6 and the M-M pairs have a maximum of 
(cos^(6'p)) ~ 0.7. The corresponding values for the correlation angle alignment function are 
~ 25% lower. 

These results raise the possibihty that clusters of all masses form with similar large 

alignments. We divide clusters into six mass bins, and for each mass bin identify the ap- 
proximate range of redshifts in which these clusters formed. We take this range to be the 
range of redshifts over which " ~ 10, the approximate maximum value of for each 

mass bin. Here, n is the number density of clusters of the given mass in each redshift shell. 
Figure 12 shows the resulting correlation angle alignment function for each mass bin. The 
range of masses and redshifts for each bin are shown in the figure. We find that clusters of 
all masses show strong alignment shortly after formation, with an alignment function that 
is nearly independent of mass. 

If cluster alignment is a consequence of cluster formation along filamentary superstruc- 
tures, then the alignment of two clusters is an indication of connecting structures. For each 
pair of clusters with separation r < 50/i~^Mpc we compute N^, the number of clusters 
inside a cyfinder of radius 1.5 Mpc along the fine connecting their two centers, and also 
the volume Vij inside this cylinder. The number density of clusters in these "filaments", 
Ufii = T^NijfEVij, can be compared to the mean number density of clusters, n, in any red- 
shift range. The resulting number overdensity in the filaments, nfu/n, is plotted as a function 
of the mean pointing angle cos^(^p) (the average of the two pointing angles of the pair) in 
Figure 13, for different redshifts. We find that overdensity in the connecting filament rises by 
a factor ~ 5 at all redshifts as the afignment increases from (cos^(6'p)) = to {cos^{6p)) = 1, 
indicating that aligned cluster pairs at separations r < 50 h"^ Mpc are more fikely to be 
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connected by filamentary superstructures than unaligned pairs at the same separation. The 
high overdensity (nfn/n) at any angle indicates that any two clusters within 50 Mpc are 
likely to be part of some supercluster. The overall amplitude of the overdensity is closely 
related to the cluster correlation function, and rises with increasing redshift as expected 
(e.g., Bahcall et al. 2003). Similar results are also seen using the cluster mass overdensity, 
or using the fraction of pairs that have another cluster along the connecting filament. 

5. Discussion 

We use a large-scale cosmological simulation of the current best-fit ACDM cosmology 
to investigate the cllipticities and alignments of ~ 10^ clusters of galaxies. We study the 
dependence of cluster cllipticity and alignment on cluster mass, radius, and redshift. We 
find that clusters are elliptical, with mean cllipticities (e) > 0.3. The derived distribution of 
ellipticities at 2; ~ is consistent with current observations (West 1989; Phonis, Barrow, & 
Prenk 1991; de Theije, Katgert, & van Kampen 1995; Basilakos, Plionis, & Maddox 2000). 

Our large sample of clusters enables a detailed study of the distribution of cluster 
ellipticities as a function of various cluster properties: mass, radius, and redshift. We find 
that the mean cluster ellipticity increases with cluster mass, i.e., massive clusters are more 
elliptical than less massive clusters. We also find that cluster eUipticity decreases with radius; 
cluster cores are less elliptical than cluster outskirts. Both of these trends can be understood 
in terms of infall and mergers; these affect massive clusters and the outskirts of clusters more 
significantly than less massive clusters and cluster cores. 

We find that the mean cluster ellipticity increases monotonically with increasing redshift 
from 2; = to 2; ~ 3, consistent with previous observation and simulation at low redshifts 
(Melott, Chambers, & Miller 2001; Plionis 2002; Jing & Suto 2002; Suwa et al. 2003; Rahman 
et al. 2004). This is likely due to hierarchical cluster formation and the high rate of infall 
and merger at early times. The evolution of cluster ellipticity with redshift does not depend 
on cluster mass or radius; we find ^ ^ 0.05 for radii r > 1 Mpc. Below a radius 
of 1 h^^ Mpc, the dependence on redshift weakens as clusters become more isotropic closer 
to the cluster center. Clusters are most eUiptical at early times and become less eUiptical 
at present; this can be understood as a relaxation process in a low-density cosmological 
background, with a reduced merger rate at later times (e.g.. Floor et al. 2003). Our results 
suggest that once clusters form, a period of relaxation exists without significant new infall 
or merger, allowing clusters to become more spherical. 

We find that clusters are ahgned with one another to separations of at least ~ 30 Mpc, 
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using both the pointing angle and the correlation angle methods for measuring alignment. 
The pointing angle measure shows a stronger signal at all separations, and a la— 2a alignment 
detection to ~ 100 Mpc (compared to ~ 30 Mpc for the correlation angle measure). 
The pointing angle emphasizes the alignment of clusters with filamentary structures along 
which nearby clusters would form, and the correlation angle, while similar, emphasizes the 
direct alignment of clusters with one another; the greater strength of the pointing angle signal 
suggests that alignment is related to filamentary structures. If individual cluster alignments 
are independently perturbed, as is likely, these perturbations will dilute the correlation angle 
signal more strongly than the pointing angle signal. 

The alignment increases with decreasing cluster separation, revealing a smooth align- 
ment function. The alignment function shows a steady increase in alignment strength from 
large separations of ~ 100 h"^ Mpc with the pointing angle measure (or ~ 30 h"^ Mpc, with 
the correlation angle measure) to a separation of r ~ 3 Mpc, where maximum alignment 
is seen. At smaller separations, the alignment measure decreases as clusters begin to overlap 
and alignments and ellipticities cannot be reliably measured. Distortions due to mergers and 
difficulty distinguishing elements of two clusters at such small separations render this region 
less meaningful. The alignment is independent of the cluster radius used in the analysis. We 
find that cluster alignment is correlated with the existence of filamentary structure along the 
line connecting a pair of clusters, with the cluster number overdensity increasing by a factor 
of ~ 5 as the alignment increases from no alignment to maximum alignment. The highly 
aligned clusters are therefore more typically members of filamentary superclusters. 

We find that cluster alignment evolves with redshift. The alignment increases for all 
cluster separations with increasing redshift from z — to z 3. Objects are most aligned at 
early times; we find that (cos^(0))inax decreases from ~ 0.55 at 2; ~ 3 to ~ 0.4 at z ~ 0. Our 
results suggest that alignment is a consequence of initial structure formation, and decreases 
rapidly over a timespan Az ~ 0.5 — 1 after which alignments continue to become more 
random, but at a slower rate. This supports the picture that alignment stems from early 
large-scale structure, with cluster formation by hierarchical clustering in which material falls 
into the cluster along large-scale filamentary structures (e.g., Onuora & Thomas 2000). This 
evolution of cluster alignment with redshift provides a new test of the cosmological model, 
even if present {z — 0) ahgnments do not significantly differ between models. 

Cluster alignment increases with cluster mass. The alignment extends to larger scales 
as the paired masses increase from low to medium to high- mass (L-L to M-M to H-H pairs). 
The alignment weakens somewhat for mixed-mass pairs (at all redshifts). At high redshifts 
z >1, alignments of high-mass clusters increase faster than ahgnments of low- mass clusters. 
We find that clusters of all masses form with similar strong alignments, and that alignment 
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rapidly declines to a steady value over a timespan ~ 0.5 — 1. Most low-mass objects 
are formed early and thus decline to a steady value ((cos^(^)) ~ 0.42 — 0.44) before high- 
mass objects can form and do so. This rapid evolution may be the result of a period of 
violent relaxation during which dissipation and gravitational instability increase isotropy 
and randomize cluster alignments. The strong evolution of the alignment signal seen shortly 
after cluster formation suggests that observations of cluster alignments at the characteristic 
redshifts of their formation, for the relevant cluster masses, may discriminate more strongly 
between cosmological models than observations at z — 0. 

Cluster alignments and cllipticitics both show dependence on cluster mass and redshift. 
The distribution of ellipticities and alignments can inform our models of large-scale structure 
formation, and can provide a new and independent test of the cosmological model. Our 
results provide detailed predictions of the current cosmology for the ellipticity distribution 
and alignment function of clusters as a function of cluster mass, radius, and redshift; these 
predictions can be used for comparison with upcoming observations of large cluster samples 
to high redshifts. 
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Fig. 1. — Density contours at 11% intervals show the three dimensional vs. the two di- 
mensional (projected) ellipticity (obtained from the best-fit ellipse to the 2-D projection of 
particle positions) for all clusters with mass MpoF > 2 x 10^^ h~^ Mq (at all redshifts). The 
correlation follows e ~ e2d- 
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Fig. 2. — Cluster ellipticity as a function of cluster mass in six redshift bins, with density 
contours at 11% intervals. 
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Fig. 3. — Histograms of the ellipticity distribution at the hsted redshifts. The sohd histogram 
is for the three-dimensional ellipticities (e), and the dashed histogram is for the projected 
two-dimensional ellipticities {e2d)- The dotted line shows the mean 3-D ellipticity within 
the redshift step. Error bars shown assume Poisson errors {a = Xj^Jn). The curve in the 
upper left shows the best-fit Gaussian to the observed distribution of 2; ^ Abell cluster 
ellipticities (Plionis, Barrow, & Frenk 1991); the observed clusters are more massive than 
the threshold used here, causing the slight shift to higher ellipticities. 
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Fig. 4. — Evolution of mean cluster ellipticity as a function of redshift for different mass 
clusters (Mc/, within the FOF radius). The bold line is the 3-D mean ellipticity (e) of clusters 
of all masses (M > 2 x 10^^) (the difference between the 3-D and 2-D mean ellipticities 
is negligible). The lower lines show the relationship for low- mass (2 x 10^^ < M < 5 x 
10^^ Mq) clusters, middle lines for medium-mass (5 x 10^^ < M < 10^^ Mq) clusters, 
and upper hnes for high-mass (M > 10^^/i~^Mq) clusters. For each, the sohd line is the 
mean three-dimensional ellipticity (e) , and the dashed line is the projected two-dimensional 
ellipticity (e2d). 
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Fig. 5. — Mean cluster ellipticity as a function of redshift for different cluster radii (i?FOF, R = 
1.5, 1.0, and 0.5/i^^Mpc). Solid and dotted lines represent the 3-D and 2-D ellipticities, 
respectively. The top panel shows the results for all clusters massive enough to ensure 
-Rfof > 1-5 Mpc, the middle panel for Rfof > 1-0 Mpc, and the lower panel for 
-Rfof > 0.5 Mpc. (The shght increase of 2-D over 3-D eUipticity for i?o.5 is due to the 
secondary ellipticity becoming more dominant on small scales, where the primary ellipticity 
is small; the 2-D distribution also includes projected material outside the 3-D radius, which 
can further increase ellipticity.) 
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Fig. 6. — The 3-D cluster alignment function (i.e., alignment as a function of halo separation) 
using the correlation angle method is shown for six redshift bins. Error bars represent la 
uncertainties. The alignment measure, (cos^(^^c)), is determined within bins of ±0.1 log]^Q(r), 
where r is the cluster pair separation in h^^ Mpc. The horizontal line indicates the expected 
value for random orientations, (cos^(6'))ran = |- 
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Same as Figure 6, but for the 2-D correlation angle alignment function. 



-24- 





00 1 10 

Cluster Separation (Mpc/h) 

(Full 3-D) 



100 



Fig. 8. — Same as Figure 6, but for the 3-D alignment function using the pointing angle 
method. 
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Fig. 9. — Same as Figure 6, but for the 2-D alignment function using the pointing an- 
gle method. The horizontal line indicates the expected value for random orientations, 
(cos^(^))ran = |- The points show observed ahgnments (converted to (cos^(6'p))) from West 
(1989) (triangles) and Plionis (1994) (squares), for Abell and Schectman clusters at z 0. 
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Fig. 10. — The 3-D cluster alignment function (pointing angle method) in six redshift bins 
(as in Figure 8), for different mass cluster pairs (see § 4). From darkest to lightest, the curves 
represent H-H pairs, M-M pairs, and L-L pairs, as labeled. At high redshift, the number of 
high-mass cluster pairs is too small for reliable analysis. Mixed-mass pairs are not shown; 
they have similar, but somewhat weaker, alignment functions. 
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Fig. 11. — Same as Figure 10, but for the 2-D cluster alignment function (pointing angle 
method) . 
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Fig. 12. — The 3-D cluster alignment for clusters of different masses shortly after their 
formation. Clusters are divided into six mass bins, and for each mass bin, cluster pairs are 
examined within the characteristic redshift range of their formation (see § 4) . The mass and 
redshift ranges are shown in each panel (with M* = 10^^ M©). 
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Fig. 13. — Mean number overdensity of clusters within 1.5/i~^Mpc of the hne connecting 
cluster pairs is plotted as a function of pair alignment, for all pairs with separations r < 
50 Mpc. The alignment, (cos^(^p)), is the average of the two pointing angles for a given 
pair of clusters, taken in bins of width 0.1. Redshift increases from bottom to top, as labeled. 



